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Abstract 

We determine from Polyakov loop correlators the screening masses in the de- 
confined phase of the (3+l)d SU(3) pure gauge theory at finite temperature near 
the transition, for two different channels of angular momentum and parity. Their 
ratio is compared with that of the massive excitations with the same quantum num- 
bers in the 3d 3-state Potts model in the broken phase near the transition point 
at zero magnetic field. Moreover we study the inverse decay length of the correla- 
tion between the real parts and between the imaginary parts of the Polyakov loop 
and compare the results with expectations from perturbation theory and mean-field 
Polyakov loop models. 
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1 Introduction 



The essential features of confinement, of the deconfinement phase transition at finite 
temperature and of the deconfined phase in QCD can be studied without loss of relevant 
dynamics in the regime of infinitely massive quarks or, equivalently, in the SU(3) pure 
gauge theory. 

The (3+l)d SU(3) pure gauge theory at finite temperature undergoes a confine- 
ment /deconfinement phase transition associated with the breaking of the center of the 
gauge group, Z(3) pfl [2|. This transition is weakly first order [3], HJ 13 EJ [7] and the order 
parameter is the Polyakov loop [fl [2] . 

In general, the deconfinement transition in the (d+ l)-dimensional SU(N) pure gauge 
theory at finite temperature can be put in relation with the order /disorder [8J phase 
transition of the d-dimensional Z(N) spin model through the Svetitsky-Yaffe conjecture [9]. 
According to this conjecture, when the transition is second order the gauge theory and 
the spin model belong to the same universality class and therefore share critical indices, 
amplitude ratios and correlation functions at criticality. There are so many numerical 
evidences in favor of this conjecture - see, for instance, Refs. [TOJ, [TTJ [T21 dSl EEH ESI ES] 
and, for a review, Ref. [T7] - that there is no residual doubt about its validity 

In the case of first order phase transitions, as for the (3+l)d SU(3) pure gauge theory, 
the Svetitsky-Yaffe conjecture is not applicable in strict sense, since the correlation length 
keeps finite at the transition point and a universal behavior, independent of the details 
of the microscopic interactions, cannot be established. Therefore there are no a priori 
reasons why the critical dynamics of (3+l)d SU(3) and of its would-be counterpart spin 
model, 3d Z(3) or 3-state Potts model [TH [19], should look similar. It turns out, however, 
that for both these theories the phase transition is weakly first order, which means that 
the correlation length at the transition point, although finite, becomes much larger than 
the lattice spacing. This leads to the expectation that the two theories may have some 
long distance aspects in common and that some specific issues of SU(3) in the transition 
region can be studied taking the 3d 3-state Potts model as an effective model. 

Since there are no critical indices to compare, a test of this possibility could be based 
on the comparison of the spectrum of massive excitations of the 3d 3-state Potts model 
in the broken phase near the transition at zero magnetic field with the spectrum of the 
inverse decay lengths of Polyakov loop correlators in the (3+l)d SU(3) gauge theory in 
the deconfined phase near the transition. If universality would apply in strict sense, these 
spectra should exhibit the same pattern, as suggested by several numerical determinations 
in the 3d Ising class [201 I2H M> [23] . 

In this paper we determine the low-lying masses of the spectrum of the (3+l)<i SU(3) 
pure gauge theory at finite temperature in the deconfined phase near the transition in 
two different sectors of parity and orbital angular momentum, + and 2 + , and compare 
their ratio with that of the corresponding massive excitations in the phase of broken Z(3) 
symmetry of the 3d 3-state Potts model, determined in Ref. [24] . 

In fact, we extend our numerical analysis to temperatures far away from the transition 
temperature T t in order to look for possible "scaling" of the fundamental masses with 
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temperature. Moreover, we consider also the screening masses resulting from correlators 
of the real parts and of the imaginary parts of the Polyakov loop. These determinations 
can represent useful benchmarks for effective models of the high-temperature phase of 
SU(3), such as those based on mean- field theories of the Polyakov loop, suggested by 
R. Pisarski [25] . 

The paper is organized as follows: in Section[2]we describe the methods used to extract 
the screening masses through Polyakov loop correlators in (3+1) d SU(3); in Section [3] we 
present our numerical results and compare them with expectations from the 3d 3-state 
Potts model and from Polyakov loop effective models; in the last Section we draw our 
conclusions. 

2 Screening masses from Polyakov loop correlators 

In this Section we define the general procedure for extracting screening masses in (3+1) d 
SU(3) at finite temperature from correlators of the Polyakov loop, 

P(x, y, z) = Tr JJ U 4 (x, y, z, n t a) , (1) 

n t =l 

where N t is the number of lattice sizes in the temporal direction and a is the lattice 
spacing. 

The point-point connected correlation function is defined as 

r<o(r) = (P}Pio)-(P})(Pio) , (2) 

where % and %q are the indices of sites and r is the distance between them. The large-r 
behavior of the point-point connected correlation function is determined by the correlation 
length of the theory, £ , or, equivalently, by its inverse, the fundamental mass. In order to 
extract the fundamental mass it is convenient to define the wall-wall connected correlation 
function, since numerical data in this case can be directly compared with exponentials in 
r, without any power prefactor. 

The + -channel connected wall-wall correlator in the ^-direction is defined as 

G(\z 1 -z 2 \) = Re(P(z 1 )P(z 2 )i}-\(P}\\ (3) 

where 

, N X Ny 

^ = ww Yl p K°' n v a ' z ) ' ( 4 ) 

represents the Polyakov loop averaged over the xy-plane at a given z. Here and in the 
following, Ni (i = x, y, z) is the number of lattice sites in the i-direction. The wall average 
implies the projection at zero momentum in the xy-plane. 
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The correlation function ([3]) takes contribution from all the screening masses in the 
+ channel. In fact, the general large distance behavior for the function G(\z\ — z 2 \), in 
an infinite lattice, is: 



G(\ Zl -z 2 \) = }^a n e- m ^-^ , (5) 

n 

where m is the fundamental mass, while m 1; m 2 , ... are higher masses with the same 
angular momentum and parity (0 + ) quantum numbers of the fundamental mass. On a 
periodic lattice the above equation must be modified by the inclusion of the so called 
"echo" term: 

G(\ Zl -z 2 \) = J2 a n[e~ mnlzi ' Z21 +e - m " {N *-^- z ^ . (6) 

n 

The fundamental mass in a definite channel can be extracted from wall-wall correlators 
by looking for a plateau of the effective mass, 

m eS (z) = In v (7) 



at large distances. For the 2 + -channel, we used the variational method [261 123 > which 
consists in defining several wall-averaged operators with 2 + quantum numbers and in 
building the matrix of cross-correlations between them. The eigenvalues of this matrix 
are single exponentials of the effective masses in the given channel. The variational 
method is usually adopted to determine massive excitations above the fundamental one 
in the given channel. In our case, however, it turned out that this method improved the 
determination of the fundamental mass in the 2 + -channel, but did not lead to determine 
higher excitations. Our choice of wall-averaged operators in the 2 + -channel is inspired by 
Ref. [23 and reads 



N X Ny 

P n {z) = — — P(n x a, n y a, z) P{n x a + na, n y a, z) — P(n x a, n y a + na, z) 

X y n x =l n y =l 

In most cases we have taken 8 operators, corresponding to different values of n, with the 
largest n almost reaching the spatial lattice size A^. We have determined the + and the 
2 + fundamental masses in lattice units, rho+ and 7712+ , over a wide interval of temperatures 
above the transition temperature T t of (3+l)d SU(3) and have studied 

- if in some region, lying strictly above T t , a scaling law can be found for the funda- 
mental correlation length £0 = l/m +; 

- if m 2 + has the same scaling behavior as rh +; 

- how the ratio m 2 + / m + compares with the ratio of the corresponding excitations 
in the broken phase of 3d 3-state Potts model near the transition, m 2 +/m + = 2.43(10), 
found in Ref. [21] on a 48 3 lattice. 

We consider also correlators of the (wall-averaged) real and imaginary parts of the 
Polyakov loop, defined as 

G R (\ Zl -z 2 \) = (ReP(^ 1 )ReP(^))-(ReP(z 1 ))(ReP(^)), (9) 
Gj(|zi-z 2 |) = (ImP(^)ImP(z 2 )) . (10) 
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The corresponding screening masses, rh^ and rhi, can be extracted in the same way as 
for the + mass. 

We have studied the ratio mi/rriR over a wide interval of temperatures above the 
transition temperature T t of (3+l)d SU(3) and seen how it compares with the prediction 
from high-temperature perturbation theory, according to which it should be equal to 
3/2 [3"0"II3"T] . and with the prediction from the mean- field Polyakov loop model of Ref. |29j . 
according to which it should be equal to 3 in the transition region. The interplay between 
the two regimes should delimit the region where mean-field Polyakov loop models should 
be effective. 

3 Numerical results 

We used the Wilson lattice action and generated Monte Carlo configurations by a com- 
bination of the modified Metropolis algorithm [32] with over-relaxation on SU(2) sub- 
groups |33j . The error analysis was performed by the jackknife method over bins at 
different blocking levels. We performed our simulations on a 16 3 x4 lattice, for which 
fit = 5.6908(2) [31], over an interval of fi values ranging from 5.69 to 9.0. 

Screening masses are determined from the plateau of m e g(z) as a function of the wall 
separation z. In each case, the plateau mass is taken as the effective mass (with its error) 
belonging to the plateau and having the minimal uncertainty. We define plateau the largest 
set of consecutive data points, consistent with each other within la. This procedure is 
more conservative than identifying the plateau mass and its error as the results of a fit 
with a constant on the effective masses m e ^{z), for large enough z. To illustrate our 
definition of plateau mass, let us apply it to Figs. Mjo) and|4](a), which give respectively 
the m 2 + and the rhi effective masses at fi = 5.74 versus the separation z between walls 
(the cases of m 2 + and rhi are usually the most problematic for the determination of the 
plateau mass). In Fig. E](b), according to our definition, the plateau is formed by the 
points at z = 3, 4 and 5, since they are compatible each other within la; the plateau 
mass is taken as the effective mass at z — 3, since this effective mass has the smallest 
uncertainty among the three masses belonging to the plateau. In Fig. H^a), according 
to our definition, the plateau is formed by the points at z = 4, 5 and 6, since they are 
compatible each other within lex; the plateau mass is taken as the effective mass at z — 4, 
since this effective mass has the smallest uncertainty among the three masses belonging 
to the plateau. 

Just above the critical value fit we find a large correlation length, which is not of 
physical relevance. It is instead a genuine finite size effect [8] related to tunneling between 
degenerate vacua. This effect disappears by going to larger lattice volumes or moving 
away from fi t in the deconfined phase. Indeed, by increasing fi one can see that the 
symmetric phase becomes less and less important in the Monte Carlo ensemble, up to 
disappearance; then, for large enough fi, also the structure with three-degenerate Z(3) 
sectors disappears, leaving only the real sector of Z(3). This is illustrated in Fig. [1] where 
the scatter plot of the Polyakov loop is shown for three representative fi values. 
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Tunneling can occur between the symmetric and the broken phase, and between the 
three degenerate vacua of the deconfmed phase. The former survives only in a short range 
of temperatures around @ t , and at (3 ~ 5.71 has almost completely disappeared, as we can 
see from Fig[I](b). When tunneling is active, the correlation function has the following 
expression [Sj: 

G(\ Zl - z 2 \) ~ a e- mo|21 - Z21 + b e- mtlzi - Z21 , (11) 

where m t is the inverse of the tunneling correlation length and is generally much smaller 
than the fundamental mass mo and therefore behaves as a constant additive term in 
the correlation function. In (fTTj) we have taken into account only the lowest masses in 
the spectrum and omitted the "echo" terms. The dependence on mt in the correlation 
function can be removed by extracting the effective mass by use of the combination 

, N , G(z)-G(z + 1) , s 

^*> =h G(, + l)-G(» + V <12) 

For (3 > 5.71 the only active tunneling is among the three broken minima. However, 
the separation among them is so clear (see, for instance, Fig. [2(a)), that it is possible 
to "rotate" unambiguously all the configurations to the real sector (see Fig. [21(b)) and to 
treat them on the same ground. 

We have performed simulations for several values of (3 in the deconfmed phase, up to 
/3 = 9.0, with a statistics of a few hundred thousand "measurements" at each (3 value. 
Simulations have been performed on the PC cluster "Majorana" of the INFN Group of 
Cosenza. 

As discussed in the Introduction and in Section [2J our typical observable is a screening 
mass in a given channel of angular momentum and parity, identified as the plateau value 
in the plot of the effective mass as a function of the separation z between two-dimensional 
walls. The procedure to determine the plateau value has been illustrated in Section [2j A 
typical example of the behavior of the effective mass with z is shown in Fig. [3] for the + 
and the 2 + channels and in Fig. [4] for the masses rhi and rhu extracted from Eqs. and 
ffTUl) . The deviation from the plateau value at small z can be attributed to lattice artifacts 
and to the possible effect of other states with excited masses in the same channel. Our 
determinations refer only to the fundamental masses in a given channel. The behavior 
with (3 of the fundamental masses in the + and the 2 + channels is given in Fig. [5j In 
Fig. [6] we show instead the masses m# and rhj for varying (3. A summary of all our 
results is presented in Table [TJ We observe from Tabled] that m + and rhpt are consistent 
within statistical errors, this indicating that the Polyakov loop correlation is dominated 
by the correlation between the real parts. Our results for m + at (3 = 5.70 and f3 = 5.90 
(T ~ 3T t /2) can be compared with the corresponding determinations on a 4 x 8 2 x 16 
lattice of Ref. [35], which give 0.271q;q2 and 0.64(1), respectively. The ~15% disagreement 
can be taken as a measure of the systematic effects involved in the whole procedure for the 
determination of masses. We cannot make a direct comparison with the determination 
for m + at (3 = 5.93 on a 16 3 x 4 lattice of Ref. [36], since we did not perform simulations 
at this value of (3. We observe, however, that the determination of Ref. [36j, which gives 
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0.73(5) according to our definition of plateau mass value, agrees with the results for the 
two adjacent ft values in our Table, ft = 5.90 and ft = 9.0. 

We can see that the fundamental mass in the + channel, as well as ifiR, becomes 
much smaller than 1 at ftt, as expected for a weakly first order phase transition. In the 
cases of m + an d of ifiR we have made some determinations below ft t (see Figs. [5] and [6]). 
It turns out that masses in lattice units take their minimum value just at ft t , where there 
is a "cusp" in the /3-dependence. Such a behaviour was observed also by the authors of 
Ref. [37J, whose results, when the comparison is possible, agree with ours. 



Table V. Summary of the results for the fundamental masses in the + and 2 + channels, 
for rhi and ifiR and for some mass ratios. 



ft 


m + 


rh 2 + 


m 2 +/m + 


rhi 


rh R 


mj/rriR 


stat. 


5.67 


0.2870(46) 








0.2793(65) 




750k 


5.68 


0.2293(67) 








0.2226(97) 




250k 


5.69 


0.175(10) 


0.81(10) 


4.61(86) 




0.176(14) 




770k 


5.695 


0.251(12) 


1.161(38) 


4.63(38) 


0.65(10) 


0.205(37) 


3.07(87) 


320k 


5.70 


0.3088(80) 


0.94(11) 


3.04(45) 


0.883(47) 


0.307(13) 


2.87(27) 


770k 


5.705 


0.342(17) 


1.09(12) 


3.20(50) 


1.056(23) 


0.352(23) 


2.86(15) 


750k 


5.71 


0.3615(30) 


1.31(18) 


3.63(53) 


1.092(54) 


0.386(23) 


2.98(19) 


770k 


5.715 


0.4054(51) 


1.25(30) 


3.08(77) 


1.105(40) 


0.4129(61) 


2.68(14) 


250k 


5.717 


0.4301(46) 


1.452(49) 


3.38(15) 


1.252(15) 


0.4237(43) 


2.955(65) 


320k 


5.72 


0.4319(69) 


1.500(80) 


3.47(24) 


1.164(58) 


0.4296(76) 


2.71(18) 


770k 


5.73 


0.4811(92) 


1.29(17) 


2.69(40) 


1.173(39) 


0.4762(92) 


2.46(13) 


170k 


5.74 


0.5267(52) 


1.44(15) 


2.73(31) 


1.035(98) 


0.5241(44) 


1.98(20) 


770k 


5.75 


0.5592(50) 


1.695(77) 


3.03(17) 


1.06(11) 


0.5562(49) 


1.91(22) 


770k 


5.76 


0.5857(65) 


1.849(39) 


3.16(10) 


1.227(28) 


0.5880(36) 


2.087(61) 


770k 


5.77 


0.620(14) 


1.83(13) 


2.95(28) 


1.421(26) 


0.6173(73) 


2.303(70) 


60k 


5.78 


0.650(10) 


1.79(12) 


2.76(23) 


1.27(12) 


0.635(10) 


2.00(23) 


60k 


5.81 


0.6742(99) 


1.812(77) 


2.69(15) 


1.39(10) 


0.671(10) 


2.07(18) 


190k 


5.85 


0.718(10) 


1.83(12) 


2.55(21) 


1.34(13) 


0.7042(87) 


1.90(21) 


110k 


5.89 


0.697(29) 


2.21(13) 


3.17(32) 


1.299(90) 


0.686(29) 


1.89(21) 


190k 


5.90 


0.734(17) 


2.09(18) 


2.85(31) 


1.284(98) 


0.7415(87) 


1.73(15) 


215k 


9.0 


0.671(19) 


1.30(57) 


1.93(91) 


1.10(11) 


0.670(19) 


1.64(23) 


200k 



3.1 Scaling behavior and comparison with the Potts model 

It would be interesting to find a scaling law for the fundamental mass in the + channel, 
which looks like the one which holds for a second order phase transition, with a suitable 
critical exponent. There is, however, an important caveat: while the correlation length 
in lattice units diverges at a second order critical point, it keeps finite at a first order 
transition point. Therefore, any second-order-like scaling law, when applied to the region 
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near a first order phase transition, should be taken as an effective description, which 
cannot hold too close to the transition point. 

With this spirit, we have compared our data with the scaling law 

/3i - AY . ™o+(A) / 13 x 



A2 - fit) m +((3 2 ) 

where m +(/3i) and m +(/32) are the fundamental masses in the + channel at (3i and /3 2 , 
respectively. This scaling law 0131) is an approximation of the law 



T 2 -T t J m +(T 2 ) 

valid in an interval around T t short enough that the linear approximation T\—T t oc (3\—f3t 
holds. We have considered several choices of (3% and found that for each of them there is a 
wide "window" of (3 values above fit where the scaling law ( 1T31) works, with a "dynamical" 
exponent v. Our results are summarized in Table [2j For /3i=5.72 and 5.73, which lie well 
inside this "window", the parameter v is compatible with z/=l/3, suggested in Ref. [38] 
to apply to the standard correlation function. For /3i=5.72 we have calculated also the 
X 2 /d.o.f. when v is put exactly equal to 1/3, getting x 2 /d.o.f.=0.75 in the window from 
(3 = 5.715 to (3 = 5.78. In Fig. [7] we show, for this choice of (3\, the comparison between 
data and the "scaling" function with v set equal to 1/3. 



Table 2: Summary of the fits of the mass ratios m +(/?i)/m +(/5) with the function {{f3\ — 
(3 t )/ '(f3 — f3t)) v . The second column contains the largest window of (3 values for which the 
fit has a x 2 /d.o.f. lower than 1. 



Pi 


window of (3 values 


V 


rVd.o.f. 


5.75 


5.70 - 5.78 


0.3619(90) 


0.90 


5.74 


5.70 - 5.78 


0.365(11) 


0.92 


5.73 


5.695 - 5.85 


0.329(12) 


0.77 


5.72 


5.695 - 5.78 


0.354(12) 


0.73 


5.717 


5.715 - 5.81 


0.3228(94) 


0.69 


5.715 


5.695 - 5.78 


0.358(10) 


0.86 


5.71 


5.72 - 5.78 


0.3951(76) 


0.33 


5.705 


5.705 - 5.78 


0.448(22) 


0.89 


5.70 


5.695 - 5.90 


0.3141(58) 


0.94 


5.695 


5.695 - 5.90 


0.3095(67) 


0.41 



Then, we have considered the /3-dependence of the ratio m2+/m + } shown in Fig. [HJ 
We have found that this ratio can be interpolated with a constant in the interval from (3 t 
to (3 = 5.77. This constant turned out to be 3.172(65), with a x 2 /d.o.f equal to 1.085. In 
the fit we excluded the point at /3=5.695, for which the determination of rri2+ is probably 
to be rejected (see Fig. [5]^b)). If the point at /?=5.695 is included, the constant becomes 
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3.214(64) with % 2 /d.o.f =2.21. The fact that the ratio m 2 +/m + is compatible with a 
constant in the mentioned interval suggests that rh 2 + scales similarly to m + near the 
transition. This constant turns out to be larger than the ratio between the lowest massive 
excitations in the same channels in the broken phase of the 3d 3-state Potts model, which 
was determined in Ref. [24J to be 2.43(10). 

3.2 Comparison with Polyakov loop models 

We have calculated the ratio mi/mR for (3 ranging from 5.695 up to 9.0. We observe 
from Fig. [9] that this ratio is compatible with 3/2 at the largest (3 values considered, in 
agreement with the high-temperature perturbation theory. Then, when the temperature 
is lowered towards the transition, this ratio goes up to a value compatible with 3, in 
agreement with the Polyakov loop model of Ref. [29], which contains only quadratic, cubic 
and quartic powers of the Polyakov loop, i.e. the minimum number of terms required in 
order to be compatible with a first order phase transition. The same trend has been 
observed also in Ref. [57] . 

4 Conclusions and outlook 

In this work we have studied some screening masses of the (3+l)d SU(3) pure gauge theory 
from Polyakov loop correlators in a large interval of temperatures above the deconfinement 
transition. 

In particular, we have considered the lowest masses in the + and the 2 + channels of 
angular momentum and parity and the screening masses resulting from the correlation 
between the real parts and between the imaginary parts of the Polyakov loop. 

The behavior of the ratio between the masses in the + and the 2 + channels with the 
temperature suggests that they have a common scaling above the transition temperature. 
This ratio turns to be ~30% larger than the ratio of the lowest massive excitations in the 
same channels of the 3d 3-state Potts model in the broken phase. This can be taken as 
an estimate of the level of approximation by which the Svetitsky-Yaffe conjecture, valid 
in strict sense only for continuous phase transitions, can play some role also for (3+l)<i 
SU(3) at finite temperature. 

The dependence on the temperature of the ratio between the screening masses from 
the correlation between the real parts and between the imaginary parts of the Polyakov 
loop shows a nice interplay between the high-temperature regime, where perturbation 
theory should work, and the transition regime, where mean-field effective Polyakov loop 
models could apply. 
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Figure 2: (a) Typical scatter plot of the complex order parameter P for (3 larger than 
5.705 on a 16 3 x 4 lattice. There are no states in the symmetric phase, but tunneling 
survives between the three broken minima. 

(b) Same as (a) with the tunneling between broken minima removed by the "rotation" to 
the real phase. 
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Figure 3: Effective mass in the + (a) and the 2 + (b) channel as a function of the separation 
between walls on the (x,y) plane at (5 — 5.74. 
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Figure 4: Effective mass mj (a) and rriR (b) as a function of the separation between walls 
on the (x,y) plane at (5 — 5.74. 
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Figure 5: Screening mass in the + (a) and in the 2 + (b) channel as a function of (3. 
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Figure 6: Screening mass rhj (a) and ttir (b) as a function of (3. 
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Figure 8: Ratio m\jm^ as a function of (3 in the deconfined phase. The three upper 
horizontal lines represent the constant (with its error) which fits the data (see the text 
for details); the three lower horizontal lines represent the corresponding mass ratio (with 
its error) found in the 3d 3-state Potts model 241. 
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Figure 9: Ratio mj/mji as a function of f3 in the deconflned phase. The vertical line 
corresponds to the critical (5 value. 
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